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ABSTRACT 



Recent X-ray observations of the supernova remnant IC443 interacting with molecular clouds have shown the presence 
of a new population of hard X-ray sources related to the remnant itself, which has been interpreted in terms of fast 
ejecta fragment propagating inside the dense environment. Prompted by these studies, we have obtained a deep XMM- 
Newton observation of the supernova remnant (SNR) Kes 69, which also shows signs of shock-cloud interaction. We 
report on the detection of 18 hard X-ray sources in the field of Kes 69, a significant excess of the expected galactic 
source population in the field, spatially correlated with CO emission from the cloud in the remnant environment. The 
spectra of 3 of the 18 sources can be described as hard power laws with photon index < 2 plus line emission associated 
to K-shell transitions. We discuss the two most promising scenarios for the interpretation of the sources, namely fast 
ejecta fragments (as in IC443) and cataclysmic variables. While most of the observational evidences are consistent with 
the former interpretation, we cannot rule out the latter. 

Key words. ISM: individual (Kes 69) — supernova remnants — X-rays: ISM 



1. Introduction 

The high-resolution X-ray imaging of supernova remnants 
(SNRs) interacting with dense ambient matter has shown 
that a population of previously unknown compact sources 
are sometimes present in the field of view, visible at high 
X-ray energies (above ~ 2 keV) . These are the X-ray shrap- 
nels - the ballistically moving isolated fragments of ejected 
material radiating both thermal and non-thermal X-rays 
(|Bvkov et al.l[2rjoa . [200l . 

Ejecta fragm ents ha ve been long known 
(lAschenbach et all [19951: iMivata et all 120011: 
lLaming fc Hwangp2003 ). but when they reach a dense 



clo ud, they may beco me bright and compact sources in 
and X-ray, with luminosity Lx > 
120021 120031) . Since the fragment 



IR (jBvkov et all 1200 
10 31 



erg s 



(jBvko 



lifetime in the cloud is about a few hundred years, a few of 
such isolated sources are expected to be visible in the field 
of an SNR entering a molecular cloud. The shrapncl-typc 
sources produce multicomponent X-ray spectra consisting 
of: (i) a relatively faint thermal continuum with kT 1 
keV, (ii) a power-law nonthermal continuum with photon 
index ~ 1.5, and (iii) non-thermal line emission due 
to K-shell ionization produced by the intense flux of 
electrons (accelerated at the bow-shock) that collide with 
a metal-rich ejecta fragment. 



The ballistically moving clumps are potentially an im- 
portant source of information about the explosive super- 
nova (SN) event. In fact, they presumably correspond to 
initial ejecta inhomogeneities, located in the faster outer 
layers, and their chemical composition may be indicative of 
the nucleosynthesis processes which occurred at early stages 
of SN evolution. Therefore, they provide complementary in- 
formation with respect to the rest of X-ray emitting ejecta. 
It should be noted, however, that a firm detection of an 
ejecta fragment is a challenging task, because fragment's 
lifetime is only a few hundred years, and the brighter is the 
fragment, the shorter is the lifetime. 

SNRs interacting with molecular clouds are the nat- 
ural sites to study the physics of fast fragments. One 
of the best cases is IC 443, which emits bright molec- 
ular lines of OH, CO, and H2, excited by a shock. 
IC 443 is an extended object (about 45' in size at a dis- 
tance of about 1.5 kpc) with the brightest GeV gamma- 
ray emitting SNR shell (lAbdo et all |2010D. It is also a 
source of TeV gamma-rays ([Albert et all 120071 ). Hard X- 
ray sources in IC 443 related to shocked clumps have 
been dete cted by Bocchino fc Bvkqvl (l2000l) and furthe r 
studied bvlBocchino fc Bvkovl (120031) iBykov et all (|2005l ). 
iBocchino et all (|2008l) and lBvkov et all (|2008l ). " 

Motivated by the results obtained in the SNR IC 443, 
we have started an observational campaign aimed at de- 
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tection of fast fragments in SNRs interacting with molec- 
ular clouds. The SNR G21.8-0.6 (Kes 69) is a good 
candidate where such sources may be seen. It is an 
exten ded incomplete radi o shell of about 20' in size 
(e.g., IShaver fc Goss 19701) at an estimated distance o f 
about 5.2 kpc (e.g- lTian fc Learry||2008t IZhou et al.ll2009T) . 
The Imaging Proportion al Coun t er ab oard the Einstein 
Observatory was used by ISewar d (1990) to obtain a noisy 
X-ray image of Ke s 69. U sing ROSAT PSPC observations, 
lYusef-Zadeh et ail (|2003t ) estimated an absorbing column 
of Nh= 2.4x10 22 cm" 2 . 

Spitzer observ a tions of Kes 69 by iReach et al.l (|2006l ) 
and iHewitt et al.l (|2009l ) have revealed bright molecular 
emission lines of OH, CO, and H2, excited by a shock. In 
particular, the SNR contains a filament along the bright 
southern radio shell of ~ 15' extension seen with Spitzer 
IRAC in the 4.5 band and attributed to the molec- 
ular hydrogen emission line, as well as an extended OH 
(1720 MHz) emission regi on that indicates t he presence 
of a molecular shock (see IHewitt et al.l I2009L and refer- 
ences therein). Millimeter band observations of CO and 
HCO+ lines towards Kes 69 performed by IZhou et"aT] 
(|2009t) provided strong evidence for the association between 
Kes 69 and the ~ +85 km s _1 component of molecular 
gas. 

The established interaction of Kes 69 with a molecu- 
lar cloud makes this SNR a promising target for a search of 
shrapnel X-ray sources. Since the size of the remnant (about 
20') is comparable with the field of view of XMM-Newton 
EPIC cameras, we conducted XMM-Newton observations 
of the field to study the faint population of hard compact 
X-ray sources and look for the signature of fast ejecta frag- 
ments. Here we report the detection of 18 sources in the field 
of Kes 69 , a significant excess with respect to the expected 
Galactic field sources of the same flux. We have found some 
of these sources to emit hard power law spectra with pos- 
sible K Q X-ray lines of Si, Fe, Ca, and Ti. Such spectra 
support the suggested shrapnel nature of the sources. 



2. XMM-Newton observations and data analysis 

The field of Kes 69 was observed with the XMM- 
Newton observatory (jjansen et al.l[200 Tl) on 2009 October 
8 (ObsID 060548010 1, PI: M. Miceli) fo r 60 ksec with the 
EPIC-MOS camera ([Turner et al.l 1200 ll) in the Full Frame 
Mode and with the EPIC-PN camera (|Striider et al.lfeOOll) 
in the Extended Full Frame Mode (the medium filter was 
used for all the cameras). In the course of XMM-Newton 
data reduction, we selected EPIC events with FLAG=0, 
patterns 0-4 for EPIC-PN camera, and patterns 0-12 for 
EPIC-MOS cameras. 

To filter out periods of enhanced particle background, 
the original event lists were screened by using the sigma- 
clipping algorithm of the ESAS softwar^E After the filter- 
ing, the net good time exposure of the field was reduced 
to about 54 ks. SAS v9.0.0 was used for EPIC data pro- 
cessing, and the HEASOFT 6.3 suite, including XSPEC 
v. 12. 3.1, was used for spectral fitting. 



1 http:/ /xmm.esac.esa.int/pub/xmm-esas/xmm-esas.pdf 



3. Results 

3.1. Diffuse emission 

Hard (3-10 keV) and soft (1-3 keV) X-ray maps of 
Kes 69 are shown in Figure [1] Diffuse X-ray emission is 
clearly seen on the latter map. The diffuse emission region 
is bordered by the radio e mission contours ob tained from 
the 1.4 GHz NVSS survey (jCondon et al.|[l998T) . No signifi- 
cant diffuse emission is found on the hard X-ray map, while 
several point-like sources are present. We verified that the 
brightest fluctuation of the residual diffuse emission in the 
hard X-ray map is at the level of only 3cr above the back- 
ground. The faint fluctuations generally correlates with the 
bright soft X-ray emission, so we argue that the very faint 
hard diffuse emission is the tail of the thermal emission vis- 
ible in the soft X-ray map, and we do not consider it any 
further. 

Spectral analysis of the diffuse extended emission of a 
mixed morphology SNR requires a thorough modeling of 
multiple components, both thermal (with possibly ovcrion- 
ized plasma) and nonthermal. Such a detailed analysis of 
the diffuse soft X-ray emission of Kes 69 is beyond the scope 
of this work. However, to estimate the absorption column 
toward this remnant and to compare it to the absorption 
of the compact sources, a spectrum of the extended dif- 
fuse emission region in the southern shell (shown as black 
rectangle in the middle panel of Figure [T]) has been ex- 
tracted and fitted with the absorbed thermal plasma mekal 
model. For the spectral fitting, a channel binning with at 
least 15 counts per bin was used. The fitting yielded Nh 
= (2.8±0.4)xl0 22 cm" 2 , T = (0.8±0.2) keV (errors are at 
the 99% confidence level) . The absorption column depth is 
consistent with that derived previously from the Einstein 
and ROSAT observatio ns. The higher temperatu re in the 
ROSAT fit derived by lYusef-Zadeh et al.l (|2003t) is most 
likely due to the contribution of the point sources unre- 
solved by ROSAT in the wide region analyzed. The post- 
shock temperature T ~ 0.8 keV corresponds to the esti- 
mated age of about 10,000 years, assuming the Sedov stage 
solution. The forward shock velocity can be estimated as 
~ 1,000 km s _1 . The velocities of the isolated ejecta frag- 
ments in the vicinity of the shell should be comparable to 
the shock velocity. 

3.2. Hard X-ray sources detection and spectroscopy 

The SAS task edetect_chain was applied to the 3-10 keV 
image to detect and study the population of point-like 
sources in the field of Kes 69 and to identify the sources 
likely associated with this remnant. The thermal emission 
from the SNR makes almost no contribution in the energy 
band chosen. Eighteen sources, whose properties are listed 
in Table [TJ have been detected with the likelihood above 30 
cr (marked in Figure [1]). The positions of the hard X-ray 
sources are remarkably correlated with the 1 2CO emission 
at the velocity a ssociated with the SNR (80 — 81 km s _1 , 
IZhou et al.ll2009l) . whose image and contours arc shown in 
Figure [TJ The correlation is better seen in Figure [2j which 
shows that the average brightness CO temperature at the 
location of the 18 sources is higher than the average cloud 
brightness temperature, which means that the hard X-ray 
sources tend to be preferentially found in dense regions in- 
dicated by bright CO emission (see Appendix A for details 
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Fig. 1. Upper panel: EPIC PN count map of Kes 69 ob- 
tained in the 3-10 keV band. The bin size is 4". The white 
circles indicate point-like sources detected in this image; 
the red contours denot e intensity of the 1. 4 GHz emission 
from the NVSS survey (jCondon et al.lll998h : the green con- 
tours denote the 12 CO (J = 1 - 0) intensity in the 80-81 
km s~l velocity range (see lower panel). The yellow cross 
denotes the position of the hard X-ray excess seen with 
INTEGRAL ISGRI. Middle panel: Adaptively smoothed 
(to a signal-to-noise ratio of 16) and vignctting-corrected 
EPIC count-rate image (MOS-equivalcnt counts per sec- 
ond per bin) of Kes 69 in the 1-3 keV band. The black 
rectangle is the region used for the spectrum of the diffuse 
emission. The superimposed hard X-ray sources and the 
radio contours are the same as in the upper panel. Lower 
panel: 12 CO (J = 1 - 0) intensity map in the velocity in- 
terval 80 -81 km/s (linearl y interpolated to a resolution of 
0.'24 fseelZhou et al.ll2009h. The contour levels are at 40%. 



Fig. 2. Histogram of brightness temperature values in the 
1 2CO maps in the 80-81 km s _1 velocity range (Zhou et 
al. 2009, also visible in the lower panel of Fig. [IJ which 
shows the molecular cloud interacting with the remnant 
(solid line). Overimposed, the histogram of the same tem- 
peratures at the location of the 18 hard X-ray sources 
detected in the XMM-Newton observation (dashed line). 
The r value is the Pearson correlation coefficient, indicat- 
ing a positive correlation (see appendix A for details) . The 
source peak is offset toward higher value with respect to 
the cloud mean value, indicating that X-ray sources cluster 
at brighter-than-average CO map points. 



about this method). The general correlation between the 
source position and the molecular cloud that is interacting 
with the remnant makes the sources good candidates for 
X-ray emitting ejecta fragments. 

We have extracted the EPIC spectra for all the sources 
rcbinned to achieve a signal-to-noise ratio > 4 per bin, and 
fitted simultaneously both MOS and PN spectra. The an- 
cillary response files were produced wit h the SAS ARFGEN 
task, and the EVIGWEIGHT task (|Arnaud et all 1200 ll) 
was used to correct for vignetting effects. An empty area 
in the northeastern part of the image was used as a back- 
ground region. Parameters of power law fits to the spectra 
of the point-like hard sources are listed in Table [1] 

Src 3, 5 and 18 were selected for a detailed study based 
on the possible presence of X-ray line emission lines in the 
spectrum similar t o those expec ted in fast moving ejecta 
knot according to iBvkovl (l2002f ) . The lines in Src 3 and 
5 can be interpreted as Fe line complex and as traces of 
the 1.8 keV Si line. The observed line features at 3.5 keV 
and 4.7 keV in Src 13 could originate from 40 Ca and 48 Ti, 
respectively. 

The spatial count distributions of Src 3, 5, and 13 are 
consistent with that of a point source. It should be noted 
here that with the few arcsec resolution of XMM-Newton 
the upper limit of a point-like source extension at the esti- 
mated distance of 5.2 kpc is about 10 17 cm. 

Two circular regions of 16" radius centered at Src 3 and 
Src 5 were chosen for spectral analysis. Unfortunately, both 
sources were projected onto the damaged MOS1 CCD6 
chip; therefore only PN and MOS2 events were taken into 
account. For spectral analysis of Src 13 a circular region 
of 20" radius containing 160 PN counts was chosen. A 10 
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Table 1. List of sources detected in the 3-10 keV band in the field of view of the XMM-Newton observation of the 
Kes 69 SNR 
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counts per bin PN spectrum was used to locate the posi- 
tions of line features in the spectrum. 

3.2.1. Src 5 

Src 5 was first detected bv lHands et~aH (|2004f ) in the XMM- 
Newton Galactic Plane Survey (XGPS). This source is 
brighter than Src 3 and Src 13, and its spectrum is shown 
in the upper panel of Figure [3] A simple absorbed power- 
law model provides x 2 — 63.9 (at 46 dof), but it does not 
fit the bright line-like features visible at about 1.8 keV and 
6.5 keV. The quality of the fit improves significantly by 
adding two narrow gaussians to describe the linefl We 
found that the line energies are 1.84±g;g| keV and 6.64±°;™ 
keV, respectively, the photon index is T = 0.6^q'2, and the 
improved x 2 is 47.8 (at 43 dof). Confidence contours for the 
line normalization coefficients and line energies are shown 
in the upper panels of Figure |U The line features can be 
associated with K-shell transitions of Si (at 1.85 keV) and 
Fe (at 6.64 keV) originating from a metal-rich ejecta knot. 
Assuming that the source is at the same distance of 5.2 kpc 
as Kes 69, the line luminosities arc ~ 6xl0 39 photons s _1 
for the Si line, and ~ 3xl0 39 photons s _1 for the Fe line 
(with relatively large uncertainties, as shown in Figure 21 
upper panels). Both these values and the ratio of the Si and 
Fe line luminosities are in agreement with those predicted 
for fast (~ 2700 km s^ 1 ) and small (size of 3 x 10 16 cm) 
ejecta knots having a mass of 10 -3 M Q moving i nside a 
dense (~ 10 3 cm" 3 ) molecular cloud (see Table 1 in lBvkovl 
l2002t where luminosities for a fragment are shown). 

2 Since the absorbing column Nb_=S.8~^\'^ x 10 22 cm~ 2 is con- 
sistent with that found for the whole remnant, we fixed it here 
to 4 x 10 22 cm -2 to obtain more restricted estimations of the 
other model parameters. 



We investigate a possible thermal origin of Src 5 by 
fitting its spectrum with an absorbed thermal model. 
Although the probable Si line is not well described by this 
model, the quality of the fit is still acceptable (\ 2 = 52.3 at 
46 dof). The best-fit parameter values for the thermal fit 

are fcT= 8-8±l;! keV and N B ^(7.8tl A e ) x 1022 cm " 2 - Since 
in the thermal scenario the column density is higher than 
that found for Kes 69, a lower limit for the distance to 
Src 5 is 5.2 kpc (i.e. the distance to the remnant), and its 
X-ray luminosity is >10 33 erg/s. These extreme values can 
be indicative of a very energetic coronal flare, as X-ray lu- 
minosities up to ~10 33 erg s _1 and temperatures up to 10 
keV have been observed in the m ost energetic flares of ac- 
tive stars fe.g-. lFavata e t al. 200^). Another possible source 
of such a high temperature and luminosity could be a cata- 
clysmic variable star (CV). Nevertheless, we point out that 
by comparing th e source count-rate in our observation with 
that observed bv lHands et all (|2004t) in the XGPS observa- 
tions (17.1±0.7 PN counts/ks and 19.3±3.8 PN counts/ks, 
for Src 3 and Src 5, respectively, in the 0.4-6 keV band), 
we found that the source luminosity is consistent with be- 
ing constant . Also, the hardness ratio of 0.75 obtained by 
lHands et al.l (|2004f ) is very similar to 0.70±0.05 obtained 
in our new observation. The stationarity of the X-ray lumi- 
nosity and hardness ratio both concur in making a coronal 
flaring or a CV origin of Src 5 unlikely, though still not 
excluded (see also Section [5]). 

3.2.2. Src 3 

The background-subtracted spectrum of Src 3 is shown in 
the middle panel of Figure [3l This spectrum is substan- 
tially absorbed (model Nn values are about 10 22 cm -2 - 
consistent with that of the soft diffuse emission from the 
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Fig. 4. Confidence contours of the line centroid and nor- 
malization for the Fe (upper left) and Si (lower left) lines in 
Src 5, and for the 48 Ti (upper right) and Ca (lower right) 
lines in Src 13. The contours correspond to the 68% and 
90% confidence levels. 
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Fig. 3. XMM-Newton EPIC-PN (black) and MOS2 (red) 
spectra of Src 5 (upper panel), Src 3 (middle panel), and 
Src 13 (lower panel) with best-fit models (power law plus 
two narrow gaussians) and residuals. 



3.2.3. Src 13 

A background-subtracted spectrum of Src 13 is shown in 
the lower panel of Figure [3] The spectrum contains line-like 
features at about 3.5 keV and 4.7 keV. Confidence contours 
for the line normalization coefficients and line energies are 
shown in the lower panels of Figure 2J A power law model 
with r = I.I^q q and two narrow gaussians at fixed energies 
fits the spectrum (x 2 =11.0 at 31 dof). The line luminosities 
at 5.2 kpc arc about 7.3xl0 38 photons s -1 and 1.7xl0 39 
photons s _1 for the 3.5 keV and 4.7 keV lines, respectively, 
in agreement with the range of values predicted for an ejecta 
fragment in a molecular cloud. Unlike the spectra of Src 3 
and 5, the spectrum of Src 13 does not show Fe complex 
line features, a fit with a thermal model is much worse than 
the non-thermal one. 



remnant), and contains line- like features at about 1.8 keV 
and 6.7 keV. Like in Src 5, we fixed the Ah to 4x 10 22 cm~ 2 , 
and found that a power law model with T = 1.3ig.7 an d 
two narrow gaussian^l describes the spectra significantly 
better than a simple power law model (x 2 =20.4 at 21 dof, 
and % 2 =28.5 at 22 dof, respectively). Assuming that Src 3 
is at 5.2 kpc, the Fe line luminosity is (1.6±0.9)xl0 39 pho- 
tons s _1 and the Si line luminosity is < 6x 10 39 photons s _1 . 
The best- fit values are in good agreement with those pre- 
dicted by iBvkovl (|2002D for a supersonic ejecta fragment 
propagating in a molecular cloud, the only difference from 
Src 5 being a smaller best-fit velocity (~ 1000 km s _1 ). We 
investigated a possible thermal origin of Src 3 by fitting its 
spectrum with an absorbed thermal model and found that, 
with such a model, we can only constrain the temperature 
to be above 7 keV. 



3 Because of the poorer statistics, we fixed the line energies to 
those found for Src 5. 



4. Hard X-ray INTEGRAL ISGRI Data of Kes 69 

The field of Kes 69 was a target of Galactic plane scans with 
the INTEGRAL gamma-ray observatory. The archival data 
obtained with the ISGRI camera during the scans comprise 
about 1700 ks of fully coded field of view observations in 
the years 2002 - 2010. 

A weak hard X-ray source is marginally seen in the cen- 
tral part of the remnant (see the upper panel of Figure [1]). 
The source position is close to the projected center of an 
imaginary triangle with vertices in Src 3, 5, and 13. The 
HEAVENS3 survey (|Walter et al.ll2010l) provides 3ct detec- 
tions in the 18 - 50 keV and 50 - 150 keV bands (see the 
lower right panel of Figure [3]). To estimate the hard X-ray 
flux of the sources, i ntercalibration with a bright nearby 
SNR G21.5-0.9 (e.g.. lKrivonos et al.l[2007l ) was performed 
to yield 5.7xl0~ 12 erg cm -2 s -1 in the 18 - 60 keV band, 
and 1.2 xl0~ 12 erg enr 2 s" 1 in the 60 - 120 keV band. 

The highest plasma temperatures consistent with the 
XMM-Newton data, discussed in §3, are too low to yield 

4 http:/ /www.isdc.unige.ch/heavens/ 
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Table 2. Fluxes of possible infrared counterparts to Src 3, 5, and 13. Infrared fluxes of Src o51, o31, o32, ol3 are quoted 
for Spitzer IRAC bands II, 12, 13, 14, Spitzer MIPS band Ml, and 2MASS bands J, H, K s (see Figure©. The flux units 
are 10 -13 erg cm -2 s _1 . Errors are at the 68% level. Estimates of extinction factors (i.e. the value 10 Ai>a ™ d / 2 - 5 , Ay is 
computed from the X-ray absorption) at positions of Src 3 and 5 are listed as ext5 and ext3. 



Band 


11 


12 


13 


14 


Ml 


J 


H K s 


A, /im 


3.6 


4.5 


5.8 


8.0 


24 


1.2 


1.7 2.2 


Source 
















o51 


2.2±0.1 


1.5±0.1 


<1.6 


<1.6 


<0.9 






o31 


4.5±0.1 


2.4±0.1 


<1.2 


<1.4 


<2.7 




3.8±0.4 4.5±0.3 


o32 


210±0.9 


120±0.4 91.0±0.8 


54.1±0.4 6.3±0.1 6.3±0.2 


75±1.6 140±2.7 


ol3 


5.4±0.4 


2.9±0.3 


2.0±0.3 


1.5±0.3 




2.2±0.2 6.6±0.2 <5.9 


ext5 


4.0 


3.5 


3.2 


2.8 


1.0 


575 


58 13 


ext3 


2.5 


2.4 


2.2 


2.0 


1.0 


76 


16 5.8 



ISGPJ detection above 60 keV. Thus, the detection of the 
INTEGRAL source associated with Kes 69 is evidence in 
favor of a non-thermal origin of the emission from (at least 
some of) the compact sources of this remnant. The power- 
law models for Src 3, 5, and 13 obtained from the XMM- 
Newton data would yield (1-5) x 10~ 12 erg cm~ 2 s _1 in 
the 18-60 keV band, which is compatible with the flux de- 
tected by ISGRI in the same band. The extrapolation of 
the power-law model of the three sources to the 60-120 keV 
band would yield a flux higher than the flux observed by 
ISGRI in that band. Therefore, if the sources were indeed 
the X-ray counterpart of the ISGRI emission, a spectral 
break would be required around 60 keV. 



5. IR imaging and photometry of Kes 69 

The field of Kes 69 was the target of Spitzer MIPS ob- 
servations rl5602176, rl5584768, rl5621120, rl5634432, 
and rl5627008, performed on 2005 Sep 30 (PI: S. Carey) 
and Spitzer IRAC observations rl2103680, rl2107264, and 
rl2110592, performed on 2004 Oct 07 (PI: E. Churchwell). 

The standard MOPEX 18.4.9 software (jMakovoz et al.l 
[2001 was used to construct mosaic images and detect point 
sources from the archival BCD-level data pre-processed by 
the S18.7.0 (IRAC) and S18.12.0 (MIPS) pipelines. The 
net exposure of the mosaic maps was equal to 2-5 frames 
of 1.2 s (IRAC) and 10-18 frames of 2.62 s (MIPS). 

Two point-like sources were found in the immediate 
vicinity of Src 3 (see the upper panels of Figure [5]) at 
a = 18 h 33 m 12 i !0, 6 = -10°13'49" (source o31, 6.2" from 
the X-ray source) and at a = 18 h 33 m 12!6, 5 = -10°13'53" 
(source o32, at 4.7"). Two point-like sources were found 
in the immediate vicinity of Src 5 (see the lower panel of 
Figure© at a = 18 h 32 m 32?2, 6 = -10°11'45" (source o51, 
at < 1" from the X-ray source) and at a = 18 h 32 m 36"5, S = 
-10°11'46" (source o52, at 4.3"). All these positions are 
in the J2000 reference frame. Another point-like source 
is clearly seen at about 2" from Src 13 (see Figure © at 
a = 18 h 33 m 07!2, S = -10°00'48" (source ol3). 

The source o5 2, listed in the USNO-B1.0 catalogue 
(jMonet et al.ll2003D . has B = 17.7, R = 14.1, which means 
that it is most likely a foreground star. The other four 
sources can be treated as the possible infrared counterparts 
of Src 3,5, and 13. Medium IR fluxes of these sources mea- 
sured w ith MOPEX and near IR fluxes listed in the 2MASS 
survey (jSkrutskie et al.l l2006t) are summarized in Table [H 
along with extinction factors. 



IRAC 3.6 urn 



MIPS 24 urn 




MIPS 24 urn 

2 18:32:36.C 



Fig. 5. Environments of Src 3 (top panels) and Src 5 (bot- 
tom panels) as seen by Spitzer IRAC camera at 3.6 and 
4.5 /im, and MIPS at 24 /im. The sources Src 3 and Src 5 
detected by XMM-Newton are shown as solid red circles, 
while the IR sources are shown as dashed green circles, both 
with labels with names referred to in the text. 
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Fig. 6. Close environments of Src 13 as seen with Spitzer 
IRAC at 3.6, 4.5, 5.8 and 8.0 /im. The red circle marks the 
position of the XMM-Newton detected source. 
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As no specific data can be found in the available catalogs 
to estimate the distances and velocities or somehow indicate 
the nature of the IR sources that are possibly connected 
with Src 3, 5, and 13, it is useful to compare the IR colors 
of Src o51, o31, o32, and ol3 with those typical for various 
classes of objects, paying a particular attention to ejecta 
fragment candidates and CVs, which have been pointed out 
as plausible interpretations of Src 3, 5 and 13. 

The Fn/Fi 3 and Fi 2 /Ft 4 flux ratios of Src o51, o31, o32 
and o l3 are larger than those in Figure 2 of IReach et al.l 
(|2006fh which contains examples of IRAC flux ratios of 
shock-excited zones in SNRs. The mn-mi2 and mi2-mi3 
colors of Src o51, o31, o32, o !3 are well within the 
range of IRAC colo rs reported bv lStolovv et al.l (|2006h and 
iGezari et al.l ((200(1 for faint X -ray sources in the Galactic 
Center region, which are usually interpreted as accreting bi- 
naries. The IR spec tral energy distribu tion of a CV source 
AE Aqr reported bv lDubus et al.l (|2007l ) is similar to that of 
Src o32, but 24/im flux of Src o32 rescaled to the distance of 
AE Aqr (102 pc) appears much lower. Thus, the assumption 
that the potential IR counterparts of Src 3 and 5 are CVs 
seems unlikely, although there is also no firm evidence for 
detection of an ejecta fragment in the IR bands. However, 
it should be noted that the ~10 34 erg s _1 IR luminosities of 
the possible counterparts of Src 3, 5, and 13 are within the 
range predicted for IR emission of ej ecta fragments inte r- 
acting with molecular clouds (see e.g. lBvkov et al.ll2008|) . 



6. Discussion 

6.1. Interpretation as ejecta fragments 

The apparent correlation of the location of the hard X-ray 
sources detected in the field of Kes 69 with the CO emission 
of molecular material associated with the remnant together 
with consistent values of absorbing columns obtained from 
X-ray spectra of the sources naturally leads to an SNR- 
rclatcd interpretation. Moreover, there is an apparent ex- 
cess of the sources with 2-10 keV fluxes above (2-3) x 10~ 14 
erg cm~ 2 s _1 within the region of about 15'x20' size where 
the SNR is interacting with the molecular cloud (shown in 
Figured]). The log(TV) — log(5) statistics of X-ray sources 
in the XMM- Newton Galactic plane survey bv lMotch et al.1 
(|2010f ) predicts about 40 sources of 2-10 keV flux above 
3xl0 -14 erg cm -2 s _1 per square degree. Therefore, only 
about 3-4 galactic field sources are expected in the region 
on the statistical ground. The X-ray luminosities of the 
sources are Lx > 10 31 erg s _1 if they are located at the 
estimated distance of Kes 69. 

If the source are related to Kes 69 , then isolated "shrap- 
nels" or "knots", associated with fast moving ejecta frag- 
ments, form an interesting explanation for a population of 
hard X-ray sources in the field of SNRs. A massive individ- 
ual fragment moving supcrsonically through a molecular 
cloud would have a luminosity Lx > 10 31 erg s _1 in the 
1-10 keV band and would be observable with XMM-Newton 
and Chandra from a few kpc distance. Such fragments have 
indee d been detected in hard X-rays inside the IC 443 
SNR (jBocchino et all 120081 iBvkov et al.l 120081 and refer- 
ences therein), and the X-ray sources i n Kes 69, which i s 
also interacting with molecular material (IZhou et al.ll2009h , 
can also be interpreted in the same way. 

The model of fast super nova fragmen ts predicts two X- 
ray emission components (Bykov 2003). The first one is 



thermal X-ray emission from a hot shocked ambient gas 
behind the fragment bow shock, with the spectrum of an 
optically thin thermal shocked plasma of an ISM cloud 
abundance. The second emission component is nonthermal. 
Interaction of the fast electrons (accelerated at the frag- 
ment bow-shock) with the fragment body produces both 
hard continuum and X-ray and IR line emission, in cluding 
the K -shell lines of the Si and Fe group elements (IBvkovl 
I2002M 

As mentioned in Sect. 3.2, the upper limit of a point- 
like source extension at the estimated distance of 5.2 kpc is 
about 10 17 cm. Individual compact fragments of scale sizes 
below ^10 17 cm would be optically thick to the K-shell 
X-ray line absorption if the metal content in the source 
exceeds about 10 -4 M Q . Indeed, for Li-like to F-like ions 
the probability of auto- ionization by a K a photon (i.e., the 
resonant Auger destruction) is si gnificantly lar ger than that 
of photon re-emission (sec, e.g., Liedahl 2005), resulting in 
the true line photon absorption. The line optical depth tx 
with respect to the resonant absorption of a K a line of 
energy Ex (in keV) for an element X (of atomic weight A) 
can be estimated as 

rx * 3 A.^.^-Nx.it [Te + A ■ (2.3 x 10 3 /? CT J 2 ] E x \ (1) 

where is the relative abundance of the ion i of the 
element X, Nx,i7 is the column density of the element X 
in the source along the line of sight in units of 10 17 enr 2 , 
and fi u is the relative line oscillator strength of the K Q line. 
The characteristic line width is assumed to be equal to the 
Doppler full width at half maximum, in which both the 
thermal ion velocity for temperature T — 10 6 • Tq K and 
the micro-tu rbulent velocity (0 a = o~ v /c) are accounted 
for (see, e.g.. iNicastro et allll999ji . 

This implies, in particular, that the line fluxes of 
such fragments are expected to be about 10 31 erg s _1 
for most of the X-ray lines, and that, for example, 
40 Ca and 48 Ti line fluxes are expected to be compara- 
ble even though the 40 Ca abundance typically exceeds 
that of 48 Ti in supernova ejecta (line saturation effect). 
The maximal model yields of 4 8 Ti in core-collapsed SNe 
range from abo ut 4xl0~ 5 (IWooslev fe Weaverl Il995t) 



to 2x10 4 Mq ([Thielemann et al" 1996ft . and are subject 



to rather substantial model uncertainties. Although core- 
collapsed SNe are natura lly associated with molecular cloud 
environments (see e.g. I Jiang et "all I2010T) and reference 
therein, still SN la type origin can not be excluded for 
Kes 69, since no compact remnant has been found there 
yet. Accor ding to the models of SN la nucleosynthesis de- 
veloped bv lMaeda et"aTl (|2010l ). the yield of 48 Ti can exceed 
10 -3 Mq in SN la ejecta. Therefore, one may conclude that 
some ejecta fragments can show line emission of 48 Ti at the 



As for the radio emission, the electrons accelerated at the 
fragment bow shock in most cases are in the MeV regime, well 
below GeV. So synchrotron radio emission is not expected to 
be observable for small fragments. However, an HII region may 
be created due to gas ionization by a powe rful soft X-ray - far 
UV emission of the fragment bow shock. IBvkov et al.l (2008) 
made some simple modeling of the radio continuum from such 
an HII region that is expected to be produced by a fast ejecta 
fragment moving through dense ambient matter. The scale size 
of the radio emitting HII region is below arcmin if the distance 
to Kes 69 is about 5 kpc and its flux at 1.4 GHz is likely about 
30 mjy (scaling from the IC443 case). 
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level estimated at 68% significance in Src 13 without any 
signs of 44 Ti, which is expected to decay at the estimated 
age of Kes 69. In conclusion, both the location, the X-ray 
spectra and the infrared counterpart of Src 3, 5, and 13 
are consistent with the expectations for supernova ejecta 
fragments interacting with a dense ambient medium. 

6.2. Alternative explanations 

The X-ray spectra of Src 3,5 and 13 may also be fitted with 
a very hot thermal component. We have therefore explored 
other explanations for the detected sample of hard X-ray 
sources. Accreting dwarfs and stars wit h active coronae can 
be sources of hard X-ray emission (e.g. Fabian et al.lll976t 
ICordova et al.l Il98lt iPatterson fe Raymond! 1 19851) and are 
therefore an interesting possibility. The X-ray luminos- 
ity of s ome CVs can reac h up to Lx ~ 10 34 erg s _1 
(see e.g. IStacev et al.l 1201 it ), but the majority of CVs are 
much fainter. Integrated emission of CVs, especially of in- 
termediate polars, is considered to be the dom inating com 



pone nt in both the Galactic Center (see e. g. iMuno et al 
2009 ) and Galactic Ridge em ission (see e.g. ISazonov et al 
20061 iRevnivtsev et al.l 12009^1 . Ind ividual spectra of ma; 



netic CVs, obtained with ASCA bv lEzuka fe Ishidal (|199£ 
show the Fe K a emission lines, both the fluorescent emis- 
sion at ~ 6.4 keV and from highly ionized Fe ions. The spec- 
tra can be attributed to the hot postshock plasma emission 
behind the accretion shock standing above the white dwarf 
surface. X-ray reflection from the white dwarf surface can 
account for the observed fluorescent Fe K a line. Hard X- 
rays in the 17-60 keV band have been detected so far by 
INTE GRAL IBIS (see iLebrun et ail 120031: lUbertini et al.l 



120031) from 37 CVs (see, e.g.. iKrivonos et al.ll2010D . Thus" 



CVs could be alternative candidate sources for Src 3 and 5, 
but not for Src 13, due to the putative presence of Ca and 
48 Ti lines along with the absence of Fe lines in its spectrum. 
However, while the ejecta shrapnel model can account for 
the emission in the 60 - 120 keV band with the flux of 
about 1.2 xl0~ 12 erg cm" 2 s" 1 detected by INTEGRAL 
ISGRI from Kes 69 field, the CVs most likely can not. Since 
the ejecta shrapnels are expected to be extended, dedicated 
Chandra observations can constrain the spatial extensions 
of Src 3, 5, and 13 to distinguish between the two interpre- 
tations. 

7. Summary and conclusions 

We have performed a study of the supernova remnant 
Kes 69 to search for compact hard X-ray sources. Such 
studies, carried out in remnants which are interacting with 
molecular clouds, may reveal a population of fast ejecta 
fragments, which may provide useful hints for the rem- 
nant progenitors. We have detected 18 sources in the 3-10 
keV band, about 3 times the expected number of galactic 
sources in this area of the sky. Moreover, the position of the 
source are not random, but correlate with the CO emission 
of the cloud interacting with the remnant. We have selected 
3 of the sources for further spectral analysis on the basis of 
the brightness and visual inspection of their spectra (Src 3, 
5 and 13 in Table [1]). The spectra can be fitted with a non- 
thermal model plus a few emission lines corresponding to 
K-shell transition of Si and Fe. We have also analyzed the 
Spitzcr IRAC and MIPS image of the field, finding possible 
IR counterparts for Src 3, 5 and 13. The X-ray spectra and 



IR flux ratio are consistent with a model of a fast ejecta 
fragment propagating inside the molecular cloud, making 
Kes 69 the second SNR for which an evidence exists for 
this new population of hard X-ray sources. 

Uncertainties in the X-ray modeling makes the alterna- 
tive explanation in terms of galactic cataclysmic variables 
also possible, but the excess with respect of log(iV) — log(5) 
counts, the location of the sources, the detection of the 
remnant by INTEGRAL/IBIS and some of the IR flux ra- 
tios make this interpretation less probable then the other. 
Further data are needed to say the final word about the 
nature of the sources. 
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Appendix A: A simple method to assess the spatial 
correlation between a set of point-like source 
and the CO diffuse emission maps 

Fig. [TJ upper panel, shows the detected hard X-ray sources 
and the 12 CO contours emission in the 80-81 km s _1 range, 
which is the range of velo city associated with the Kes 69 
remnant (jZhou et al.ll2009t ). By eye, the source seem to be 
correlated with the cloud emission, but one need to assess 
this correlation in more detail. A simple way to do this is 
described below. 

We have computed the histogram of the 80-81 km s _1 
12 CO brightness temperature map shown in Fig. [TJ lower 
panel. This histogram is plotted as a solid line in Fig. [5] and 
has a peak around 2 — 2.5 K. We have also computed the 
brightness temperature at the locations of the 18 detected 
hard X-ray source in the same map, and we also computed 
the histogram and overplottcd it as a dashed line in the Fig. 
[2] We stress that in case of a uniform distribution of the 
hard X-ray sources, their histogram should match the cloud 
histogram. If the sources tend to be located preferentially 
outside the cloud, the peak of their histogram should lie to 
the left of the peak of the cloud. Conversely, if the sources 
tend to be inside the cloud, their peak should fall at a tem- 
perature higher than the cloud temperature, which means 
that the source have a tendency to be located in regions of 
high brightness temperature. This is exactly what happens 
in case of the histograms computed using the 80 — 81 km 
s -i i2qq brightness temperature map (Fig. [5]). 

To further strengthen this argument, we have also com- 
puted the histogram using the 8-10, 50-60, 53-55 and 68-70 
k m s" 1 12 CO bright ness temperature maps, using the data 
of IZhou et aDl2009L and we show them in Fig. IA.1I These 
velocity ranges are not associated to the remnant and there- 
fore they should be in foreground. Indeed, the plots show 
that the source histogram either has the peak at the same 
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Fig. A.l. From top to bottom: as in Fig.[5]but for 8-10, 50- 
60, 53-55 and 68-70 km s _1 12 CO brightness temperature 
maps 



brightness temperature of the cloud histogram, or at lower 
brightness level, indicating little or negative correlation. 

We have also calculated the Pearson's correlation coef- 
ficient r between the CO image and the hard X-ray sources 
image (an artificial image with the pixel values of one at 
the source position and zero elsewhere) . This coefficient is 1 
if the images are identical, > if the two images are corre- 
late, if they are completely uncorrelated, and < if they 
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are anti-correlated. The correlation coefficient r is shown 
in the upper-left corner of each image. The 80 — 81 km s _1 
histogram is the only one with positive r, i.e. the only one 
indicating a possible correlation between the sources and 
the 12 CO map. 
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